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Abstract 

Atomic and molecular electric dipole moments are calculated within the mini- 
mal supersymmetric standard model. Present experiments already provide strong 
bounds on the combination of phases responsible for the dipole moments of the 
neutron and closed shell atoms. For a supersymmetry breaking scale of 100 GeV, 
these phases must be smaller than ~ 10 -2 . 
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The experimental bounds on electric dipole moments (edm's) of the neutron, 
atoms, and molecules, are reaching a level of precision that makes them useful 
probes of P and T violation that may originate at scales beyond the standard 
model. Aside from the QCD vacuum angle, the only source of T (or CP) violation 
in the standard model is the phase in the quark mass matrix. This generates edm's 
orders of magnitude beyond experimental reach |TJ]. Extensions of the standard 
model however generally have CP violating phases that produce edm's of the order 
of, or larger than, present experimental bounds. A finite 6qcd that saturates the 
bound from the neutron edm 0, would also produce atomic edm's of the same 
order. In this paper the main sources for atomic and molecular edm's in the 
minimal supersymmetric extension of the standard model will be identified. 

In order to identify the microscopic sources of P and T violation responsible 
for the edm's, it is instructive to present the analysis starting from the largest 
relevant scale, namely the atomic scale. The next step is then to proceed down 
to the nuclear scale. This will enable finally a discussion of the origin of P and T 
violation at the subnuclear level. 

At the atomic scale an atom or molecule is a composite system of electrons and 
nuclei. These constituents contribute to the edm via T and P odd electromagnetic 
moments and local non-electromagnetic interactions. Some of the electromagnetic 
moments are suppressed. In the nonrelativistic limit, the contributions from the 
electron and nuclear edm's vanish. This is known in the literature as Schiff's 
theorem ||. The electron edm can contribute though in heavy atoms where the 
electrons are relativistic j|, ||. But in closed shell atoms and molecules with paired 
electron spins, the electron edm contributes only through hyperfine interactions 
with the nucleus, and is therefore suppressed ||. Schiff's theorem does not apply 
to higher moments such as the nuclear magnetic quadrupole moment. In addition 
there is another T and P odd moment known as Schiff's moment, which is pro- 
portional to the offset of electric charge and dipole distributions of the nucleus. It 
leads to a local electromagnetic coupling of the electron to the nucleus [Q. 
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At the nuclear scale all of the contributions discussed above are contained in 
the following Hamiltonian 

H = l -d e ea^i lb eF^ + ^d N Na^i lB NF^ + C P n s>s ^(ei l5 e)(NN) 

+dy s ^(ee)(Ni l5 N) + C^ PT ^(ea, u e)(Na^i l5 N) 

+Qv,Pv(ei,e)(N S> ^N) + Q e P n vy (e % l5 e){N^N) 

+ C n s n PS ^=(NN)(M l5 N) + Q n v y v (N 7fl N)(N # l5 N) (1) 

where N = (p, n) and isospin violation is ignored. The first and second terms are 
the electron and nucleon edm's respectively. The remaining terms are local inter- 
actions among the electrons and nucleons. The nucleon edm and nucleon-nucleon 
couplings contribute to the Schiff and magnetic quadrupole moments 0. The 
electron-nucleon couplings of course contribute to the electron-nucleus coupling 
||. Calculations of the contributions of the terms in H to the edm's consid- 
ered below may be found in Ref.s [5,7-11] and are summarized in Table 1. The 
contributions from Q\? PV , Q e pvvi an d Qvpv have not been considered in these 
references although Qy" P y does contribute directly to the Schiff moment. These 
derivative operators have the same nonrelativistic limit as some of the nonderiva- 
tive operators. This allows a rough estimate of the corresponding contributions to 
the edm's. 

In this paper the microscopic origin of the operators in (1) will be calculated 
within the context of the minimal supersymmetric extension of the standard model 
| 12|1 . In addition to the KM like phases, this model has other CP violating phases 



in the superpotential and soft supersymmetry breaking terms. The superpotential 
contributes a phase from the term \ix\e l ^^H\H2- The phases in the soft breaking 
terms arise from the gaugino masses, — \\ma\e l ^~ a \a\a + h.c, and the mass param- 
eters coupling left and right handed squarks (sleptons), —flmf\Af\e^ A ffR + h.c, 
where m/ is the corresponding quark (lepton) mass and Af is defined by 

m/ A, = (2) 



Ap is the scalar trilinear soft breaking coupling and h f the Yukawa coupling. The 
importance of these new phases is that P and T odd operators can be generated 
irrespective of generational mixing.[] 

The supersymmetric origin of the terms in the Hamiltonian (1) will now be 
considered. 

1. The electron edm is a chirally violating operator of effective dimension 6. It 
arises at the one loop level from electroweak gaugino exchange as shown in Fig. 1. 
The sum of all such diagrams contains several unknown masses and phases [|Dj] . 
In order to display the order of magnitude, only the photino contribution will be 
considered lfL4j| , 

d e = ~ e 24 7r W ^f £ sin (^ ~ fa)f( x ) + ( 3 ) 

where x = m\/m'~, f(x) is a loop function of order unity, /(l) = 1 |l4jl , and H — • 
represents the contributions from the other electroweak gauginos. Numerically, 
for rag = m 7 Eq. (3) gives |14| 

1 in 25 (100 GeV) 2 . , , 
d e ~-lx 1(T 25 ^ — sin Ae - 07 e cm + ■ • • 4 

my\A e \ 

2. The lowest dimension operators which potentially contribute to the nucleon 



edm are the light quark electric and chromoelectric dipole moments [|T^, [T1J, and 
Weinberg's 3-gluon operator |T7|, |TB 



01 = Iqa^qF^ 

^3 — ^Jabc^afiu^b L7 cp 

where G a ^ v = \(-^ U p a G P a an d T a are the SU(3)c generators. As for the electron 
edm the quark dipole moments are effectively dimension-6. The operators 0\ and 
2 are generated at the one loop level through gluino and electroweak gaugino 



exchange |16|]. The sum of all such diagrams is again a function of several masses 
1 If the hidden sector is Polonyi and the electroweak gaugino masses vanish at tree level, the 
supersymmetric phases can be rotated away. 
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and phases. To establish the order of magnitude only the gluino and photino 
exchange will be considered 

C M = 77ZT m <?KI ( ^ sin (<K - 4>~ x )^f(y) + ^%sin(0 Ag - fa)f(z) H J Ci 



24tt yi *' m XT q TA '3 y ' m% 

\ A / 



(5) 



> A K 

(6) 

where Q q is the quark charge, y = m~/m~, z = m~/m~, / is the same function 
that appears in Eq. (3), and h(l) = 1 Ci are renormalization group correc- 



tions for the evolution of the operators from the supersymmetric scale, M, to the 
hadronic scale, \i. Using the known anomalous dimensions |L9| with a s (M) = .1 
and a s (fj,) = 4tt/6 |17J gives Ci — 1-6, and (2 — 3.6. Note that the photino con- 
tribution is down by ~ aa~ l compared with the gluino. Barring any cancellation 
among the phases, and if all the mass parameters are of the same order, the gluino 
contribution will dominate the quark dipole moments. 

The 3-gluon operator is generated at the two loop level. The largest contribu- 
tions come from integrating out the chromoelectric dipole moments of quarks with 



mass m Q > A QC d 0> [H 



8n m Q 

Keeping only the charm quark as the dominant contribution gives £3 ~ .38. 

An estimate of the edm requires an evaluation of the matrix elements of these 
operators on the nucleon. There is however no systematic approximation scheme 
to reliably calculate hadronic matrix elements of this type. Here, some empirical 
rules known as "naive dimensional analysis" that keep track of factors of 4tt and 



mass scales [|T7|, will be used. This is the most uncertain aspect of the entire 
analysis. Even so, there is no physical reason to expect a substantial suppression 
or enhancement as compared to estimates given by these rules. Using these rules, 



the nucleon edm is 



d N ~ d(//) + -^C 2 (//) + -^A x C 3 (n) 



where A x ~ ^f-w is the chiral symmetry breaking scale. Numerically, keeping only 
the gluino contributions with rriq = m~ x = rrij = m, and A u = Ad = A c = A q , Eq. 
(8) then gives 

\d N \ ~ (2.2 + 1.1 + .1) (1 ° O 3/ G ^ )2 sin(0 Ag - A ) 10- 23 e cm (9) 

m I \ A q \ 

where the the terms on the right side come from 0\, O2, and O3 respectively. 
3. The nucleon-nucleon couplings are just the T odd component of the nuclear 



force. Following Weinberg's analysis of nuclear forces pl| , it is useful to think of 
these as arising from the T odd exchange of low lying mesons and direct contact 
interactions. For illustrative purposes only the exchange of the lightest meson, 
the pion, will be considered; and an estimate of the T odd pion-nucleon coupling, 
9ttNnNttN , arising from the microscopic physics will be made. The leading con- 
tribution will come from the light quark chromoelectric dipole moment. Using 
"naive dimensional analysis" the coupling is 

9ttNN - 97vNN^ 2L C 2 (fl) (10) 
47T 

where g W NN is the usual pseudoscalar pion-nucleon coupling constant. Integrating 
out the pion, the contribution to the nonderivative nucleon-nucleon coupling is 

s,ps ~ 9ttNn9ttNN 2r^~ v 11 '' 

TTlif (jr p 

Numerically, keeping only the gluino contribution with m~ x = trig, Eqs. (6), (10), 
and (11) give 

°™ ps " 7 ( r?/ii7 sin(0A ■ ^ (i2) 

The exchange of heavier mesons will also contribute to this coupling and the 
derivative nucleon-nucleon coupling. This is not expected to substantially change 
the estimates given below based solely on the nonderivative coupling from pion 
exchange. 

4. The electron-nucleon couplings in (1) arise from two classes of operators. The 
first operators involve an electromagnetic infrared enhancement and are suppressed 
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by two powers of a heavy mass. The second class of operators are local at the 
atomic scale and suppressed by four powers of a heavy mass. 

The first class of operators arise from the effective nucleon-photon couplings 
NNF^F^ and Ni^NF^F^ . The two photons couple to electrons to produce 
an effective electron-nucleon interaction. Integrating out the photons gives an 
infrared enhancement, cutoff by the electron mass. One way in which such an 
effective interaction arises is through the pion pole as shown in Fig. 2, where 
the CP violation occurs in the pion-nucleon coupling. This contributes to the 
electron-nucleon pseudoscalar-scalar coupling 

Cps,S ~ g-KeeQirNN — 777" (13) 

where g n NN is given in (10), g nee is the effective pion-electron coupling computed 



in Ref. 22 



3a 2 m e , ( A x 



- 2^/, ln UJ < 14 > 

and the approximation that the electron is on shell has been used. Note that the 
effective electron-nucleon interaction is suppressed by only two powers of a heavy 
mass (from the light quark chromoelectric dipole moment). A similar diagram with 
the CP violation in the vertex tt°F^ v F^ v ', contributes to the scalar-pseudoscalar 
coupling. The edm's in Table 1. are however somewhat less sensitive to this 
coupling. 

The nucleon-photon couplings can also arise from the following microscopic 
operators which are not suppressed by a light quark mass 

C5 = G ailiV G^F pa F pa 

These operators arise at the two loop level. The dominant contribution to O4 
comes from integrating out the electric dipole moment of a quark with mass ulq > 
Aqcd, as shown in Fig. 3. This involves an infrared enhancement, cutoff of by 
the quark mass. The coefficient is therefore suppressed by only two powers of a 



G 



heavy mass. Similar operators have been considered in Ref. |[23|| . Here 

eg s 2 QQ 
25Qir 2 mQ 

Using "naive dimensional analysis" then 



C^ s ~ -6^m e A x ln (^) (16) 

A similar discussion applies to 0§ with the chromoelectric dipole moment of heavy 
quarks. All these operators scale like (m^M) -2 . 

It should be noted that the infrared enhancement from the photon-electron 
loop comes from momenta ~ m e . Treating the resulting interactions as local is 
somewhat dubious for heavy atoms for which p e ~ Zam e . This is not expected to 
significantly alter the estimates of the edm's given below. Among the first class of 
operators, the estimate in Eq. (13) turns out to be numerically most important, 

Cf s>s ~ 2 x 10~ 8 (1 °^f7 sin ^ - ^ 

m ql \ A q\ 

The second class of local operators that lead to electron-nucleon couplings in- 
clude 

C> 6 = ei^eqq C 7 = eeqi^q 

Og = ea^eqa^i^q Og = Eij 5 eG a ^G^ 
Oio = eeG^GT di = e lfl eq # lh q 
12 = 6 0^ "iseqYq 

Because of the chiral properties, all these operators are effectively dimension 8. 
The operators Oq, O7, and 0% are generated, after Fierz reordering, by box dia- 
gramsQof the type shown in Fig. 4. Neglecting running corrections, the coefficients 
are then related by C 6 = C 7 = ~C 8 = C. A typical contribution is of order 

a m e m g \^\ 

C ~^ GF s^4 8m ^-^ (18) 
where tan/3 = v 2 /vi. With all the supersymmetric mass parameters of order M, 

this scale like GfM~ 2 . 

2 Tree level Higgs exchange is CP conserving in the minimal supersymmetric standard model 

and will not contribute to the electron-nucleon couplings. 
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The chiral suppression of a light quark mass is avoided in the operators O d 
and O w . These can arise from the box diagrams of Fig. 4 with the light quarks 
replaced by heavy quarks, Q. The heavy quarks are then integrated out as shown 
in Fig. 5. In the limit M » m Q , 

c 2a s C 6 (m Q ) 

q 3 An itlq 

C 10 = £-^^ (19) 

There are no dimension-8 operators without a light quark mass suppression that 
give a tensor-pseudotensor Lorentz structure. There is however an effective dimension- 
10 operator, Oi% = \eo ' p.vtidabcG ' apa G p h c T G^ V ', which is generated analogously to Og 
and (9io, but with at least three gluons on the heavy quark loop. With the rules of 
"naive dimensional analysis" this operator is not much suppressed compared with 
Og and Oiq. It is therefore included in the estimate below. Using "naive dimen- 
sional analysis" the contributions of C 6 through O 10 , and 13: to the nonderivative 
electron-nucleon couplings are 

„ \ F! / A \ 3 \ 

(20) 

24 \AitJ 1 " ! ' 
where again running corrections have been neglected. 

The operators O n and 12 are generated by box diagrams of the type shown 
in Fig. 6a, 




Cu ~ — e Q -sin ^ 

An H m° 
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Ci 2 ~fe 2 gj^^sin(0,-0 A J (21) 
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There are also contributions that result from a weak interaction between the elec- 
tron(quark) and the weak-edm of the quark(electron) as shown in Fig. 6b. 

C u = ±(1 - 4|Q g |sinX)(l - 4sinX) ° r ( : (/ ' ' 



y/2 e 



C d^ 

Ci 2 = ±(1 - 4|Q (? |sin 2 ^)cot^^^ (22) 

V2 e 

where ± refers to up or down type quarks, is the electron edm arising from 
wino exchange, and Ci(/x) is given above. These operators scale like either GpM~ 2 
or M -4 . Using "naive dimensional analysis" the contribution to the derivative 
electron- nucleon couplings in (1) are 

^V,PV ~ °H 

Qpv,v = C\2 (23) 

All of the local operators turn out to be numerically somewhat less important than 
the estimate in Eq. (13). 

The atomic and molecular edm's can now be estimated as functions of the 
microscopic parameters. In order to identify the most important effects at the 
atomic scale, the contributions from the electron edm, nucleon edm, nonderivative 
nucleon-nucleon coupling, and the largest electron-nucleon coupling (i.e. from 
Eq. (13)) will be retained. Since an edm is proportional to spin, the leading 
contributions in open and closed shell atoms are different. For open shell atoms 
we consider 133 Cs and 205 T1 since good experimental bounds are available |?4|. ^| . 
Putting together the results cataloged above with the results of the atomic and 
nuclear calculations from Table 1., 



iCs 



- 1.2 sin(0 Ae -^) + 2x KT 3 sin(0 A , - X ) + 1 x lO^sin^ - X ) 



100 GeV\ 2 



+ 1 x lO- 3 sin(0 A , - fa) ) — j 10- 2S e cm (24) 
d Tl ~ ( 6 sin((f) Ae -^) + 2x 10- 4 sin(^ g -^) + 2x lO- 3 sin(0 A - X ) 



N /100GeV\ 2 

- 1 x 10- 2 sin(^ g - 4>~ x ) ) { M J 10- 23 e cm (25) 



where, for simplicity of notation, all the supersymmetric mass parameters have 
been assumed to be of order M. The quantities on the right hand side arise 
respectively from the electron edm, nucleon edm, nonderivative nucleon-nucleon 
coupling, and pseudoscalar-scalar electron-nucleon coupling. In both cases, the 
electron edm gives the dominant contribution since it is enhanced in heavy open 
shell atoms. For 133 Cs, with a nuclear quadrupole moment, the nucleon-nucleon 
coupling is only a factor ~ 10 less important than the electron edm (again assuming 
all the masses and phases are the same order). The present experimental bounds 



are 



|c? Cs | < 7.2 x 10" 24 e cm 
\dn\ < 6.6 x 10~ 24 e cm 



The 205 T1 result gives a bound on the phases contributing to the electron edm of 

/100 GeVV 
V M J 



.100 GeV\ 2 
sin A „- fall — < -1 



For closed shell atoms, good experimental bounds are available for 129 Xe and 
199 Hg. Combining our results with those from Table 1., 

d Xe ^ ( 8 x lO- 3 sin(0A e - fa) - 7 x KT 2 sin(0 A , - fa x ) + 3 sin((f) Aq - fa x ) 

+ 2 x lO- 4 sin(0 A , - fa x ) ) ( 10 °^ eV ) 2 ID" 26 e cm (26) 
d Ug ^ ( - 1.2 x lO~ 2 sin(0 j4e - fa) - 3 x lO^sin^ - fa x ) + 4 sm{<p Aq - fa x ) 
+ 2 x lO- 4 sin(0 Aq - fa x ) ) (^^) 2 10- 25 e cm (27) 

The electron edm is here suppressed since the electron spins are paired. The 
leading contribution comes from the nucleon-nucleon coupling. The present ex- 
perimental bounds are [26-29], 

|dxe| < 1-4 x 10~ 26 e cm 

|d Hg | < 3 x 10" 27 e cm (95% C.L.) 



3 Unless stated otherwise, the experimental bounds given here are the sum of the reported 
measurement and experimental error. To date, all measurements are consistent with zero. 
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The 199 Hg result gives a bound on the phases contributing to the light quark 
chromoelectric dipole moment, 

, . /100 GeV\ 2 
sm(fa q - fo) ( — j <.008 

It should be noted that, due to the uncertainties in the hadronic matrix elements 
and nuclear calculations J7|, this bound is much more uncertain than that from 
the electron edm given above. 

At present, the best bound on a molecular edm is for 205 T1F. Again combining 
our results with Table 1., 



d TlF ~ ( - 8 x 10 2 sin(0 Ae - fa) - 1.5 x 10 L sm((p Aq - fa x ) + sin(0 j4q 

/100GeV\ 2 w 

I 1 1 n— 22 



- 1 x 10- 3 sin(<^ - X ) ) — j 10~ 22 e cm (28) 

The leading contribution again comes from the nucleon-nucleon coupling since the 
electron spins are paired. The current experimental bound [110, [3T 



|g?tif| < 4.6 x 10 23 e cm 

does not yet yield a substantial bound on the phases if M ~ 100 GeV. 

Because of the different structure of open and closed shell atoms, the exper- 
iments bound different microscopic CP violating parameters. Open shell atoms 
are sensitive to phases in the weak gaugino sector contributing to the electron 
edm. Closed shell atoms and molecules with paired electron spin get the leading 
contribution from nuclear effects. This provides a stringent bound on the gluino- 
squark phases contributing to the light quark chromoelectric dipole moment. For 
comparison, the neutron edm receives contributions from both the light quark 
electric and chromoelectric dipole moments The experimental bound of 
\d n \ < 8 x 10 -26 e cm, with the estimate (9), gives 



/100 GeV V 
V M ) 



This is of the same order as the bound from 199 Hg. Taken together, the experimen- 
tal bounds provide information on different combinations of CP violating phases. 
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Other sources of CP violation would give rise to different patterns of edm's [ SI] . 
For example, a nonzero Oqcd would contribute to the edm's mainly through the T 
odd pion-nucleon coupling J?], |32|. This would give a definite relation among the 
edm's of open and closed shell atoms and the neutron. 

The prospects for improving the present measurements are encouraging. The 
ultimate sensitivity of the current 199 Hg experiment is expected to reach the 
level of 3 x 10~ 28 e cm [29|. New techniques may allow further improvement 



in sensitivity for open shell atoms ||33|| . An unexplored area where the phases 
contributing to the light quark chromoelectric dipole moment could be measured 
is in light atoms. In this case the electrons are nonrelativistic and the primary 
contributions are nuclear. In summary, the experimental study of atomic edm's 
provides a promising probe for exploring physics beyond the standard model. 

We would like to thank D. Heinzen and V. Kaplunovsky for useful discussions. 
We would also like to thank N. Fortson for providing us with the unpublished 
bound on 199 Hg. This research was supported in part by the Robert A. Welch 
Foundation and NSF Grant PHY 9009850. 
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133 Cs 


205 rpj 


129 Xe 


199 Hg 


205 T1 p 


d/d e 


120 


-600 


-.8 x lO- 3 


1.2 x 10- 2 


80 


d/dN 


7 x 10" 4 


7 x 10" 5 


2 x 10~ 5 


~ io- 4 


.5 


d(e cm) /C&g 


2 x 10~ 24 


3 x 10~ 26 


5 x 10~ 26 


6 x 10~ 25 


2 x 10~ 22 


d(e cm)/C? St8 


7 x 10~ 19 


-5 x 10~ 18 


9 x 10~ 23 


1 x 10~ 21 


-6 x 10~ 18 


d(e cm) /ClppT 


9 x 10~ 21 


5 x 10~ 21 


5 x 10~ 21 


-6 x IO" 20 


1 x 10~ 16 


d(e cm)/C|> 5 






-1 x 10~ 23 


3 x 10~ 23 


-4 x 10~ 19 



Table 1. Contributions to electric dipole moments from the Hamiltonian (1). 



1G 



Figure Captions 

Fig. 1 A typical contribution to d e from photino exchange. 

Fig. 2 The contribution to an effective electron-nucleon coupling from pion ex- 
change with the T odd pion-nucleon coupling. 

Fig. 3 The contribution to the operator G aiJLV G^ v ' F p(7 F p(T from the electric dipole 
moment of a heavy quark. Other diagrams related by gauge invariance are 
not shown. 

Fig. 4 A typical contribution to the electron-quark operators. 

Fig. 5 The contribution to the electron-gluon operators from the electron-heavy 
quark operators. Other diagrams related by gauge invariance are not shown. 

Fig. 6 Typical contributions to the operator e d p ^eq^q; (6a) from a box dia- 
gram, (6b) from Z-exchange and the electron electroweak dipole moment. 



17 



